Fluid/Gravity Correspondence, Local Wald Entropy Current 

and Gravitational Anomaly 

Shira Chapman|*| Yasha Neimanjll and Yaron 
Raymond and Beverly Sackler School of Physics and Astronomy, 
Tel-Aviv University, Tel-Aviv 69978, Israel 
(Dated: February 14, 2012) 

Abstract 

We propose, in the framework of the fluid/gravity correspondence, a definition for a local hori- 
zon entropy current for higher-curvature gravitational theories. The current is weh-defined to first 
order in fluid gradients for general gravity actions with an algebraic dependence on the Riemann 
tensor. As a detailed example, we consider five-dimensional Einstein-Maxwell theory with a mixed 
O ; gauge-gravitational Chern-Simons term. In this theory, we construct the proposed entropy current 

. on a charged black-brane background, and show that it has a non-negative divergence. Moreover, 

a complete correspondence between the charged black-brane horizon's dynamics and the hydrody- 
namics of an anomalous four-dimensional field theory is established. Our proposed entropy current 
is then found to coincide with the entropy current of the anomalous field theory fluid. 
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I. INTRODUCTION 



The fluid/gravity correspondence relates field theory hydrodynamics to black hole (brane) 
dynamics. In the hydrodynamic regime the system is in approximate thermal equilibrium in 
each local region, while the parameters of the equilibrium vary slowly between one region and 
the next. In a particular setup the fluid/gravity correspondence follows from the AdS/CFT 
correspondence (for a review see [ij]), when considering the hydrodynamic regime of the 
field theory j^. The essential ingredients needed to relate fluids to gravity are the existence 
of a horizon in the gravitational background that is related to a thermal equilibrium state 



in the field theory, and a derivative expansion around it 



Ml. 



Thus, one can define the 



relation between fiuids and gravity on more general backgrounds, e.g. the Rindler geometry 
G-sl. One can also think about the fiuid/gravity correspondence as an extension of black 
hole thermodynamics, where charges are upgraded into local currents, and the black hole 
entropy jo], [lO| into a local entropy current. 

The local entropy current in Einstein gravity is a vector density intrinsic to the horizon 



It is directed along its generating light rays, and its fiux through a spatial slice of 
the horizon gives the area of the slice. In this work we will propose, in the framework 
of the fiuid/gravity correspondence, a definition for a local horizon entropy current for 
higher- curvature gravitational theories. This should be a vector density on the horizon 



hypersurface, such that for stationary solutions its fiux equals the total Wald entropy |12| . 
while for non- stationary solutions it satisfies an increase law (9^s^ > 0. Our definition 
of the entropy current is potentially ambiguous. However, for actions with an algebraic 
dependence on the curvature tensor, we will show that current is well-defined to first order 
in fiuid gradients. 

An interesting application is to a five- dimensional Einstein-Maxwell theory with a gauge- 



gravitational Chern-Simons term. For comp. 



eteness, we also include a pure gauge Chern- 



Simons term which was already studied in [13]. In this setup, we will construct the proposed 
entropy current to first order in gradients on a charged black-brane background, and find 
that its divergence is positive. 

This particular 5d theory has attracted much interest recently, because it implements 



chiral 
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15| and mixed chiral-gravitational 
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18| anomalies in the dual four-dimensional 



field theory. With these anomalies, it has been shown that new, non-dissipative, transport 



coefficients arise 



19|, |20|], which are responsible for the generation of current along magnetic 



fields and vortices in the fluid. A particular manifestation of these transport coefficients 
(the chiral rna^netic effect 21]) was suggested to exhibit charge separation [22| and spin 
distribution [23] in heavy ion collisions 24 1. 

To establish fully the fluid/gravity duality between the anomalous QFT hydrodynamics 
and Einstein-Maxwell theory with a gauge- gravitational Chern-Simons term, we will analyze 
the horizon dynamics along the lines of 13|], and obtain the complete fluid equations. We 
then find that the entropy current as defined in the hydrodynamics coincides with our 
proposed Wald entropy current. 

The paper is organized as follows. In section HIl we review essential ingredients for the 
discussion, including null horizons, the Bekenstein-Hawking entropy current, hydrodynamics 
and Wald entropy. In section llllt we define our Wald entropy current and study its proper- 
ties. In section llVt we consider the Wald entropy current of the 5d Einstein-Maxwell theory 
with a gravitational Chern-Simons term. In section |Vl we establish the duality between 
the horizon dynamics in this theory and the hydrodynamics of an anomalous QFT. In the 
appendices, we provide the detailed calculations behind the results of sections [TVWl 



II. PRELIMINARIES 
A. Null horizons 

We will be interested in event horizons, or more generally in null hyper surfaces. In the 
following, we briefly review some of the required ingredients for our studies. The general 
arguments are independent of the dimension, but for concreteness we will refer to a four- 
dimensional horizon embedded in a five-dimensional bulk spacetime. This is also the relevant 
case for the AdS5/CFT4 correspondence, and for the particular model in sections HVUVl We 
denote bulk coordinates by a;^, with capital Latin indices. Where relevant, we work in a 
coordinate system adapted to the horizon, i.e. we decompose = {r,x^), with the horizon 
situated at r = 0. Small Greek indices denote tensors intrinsic to the horizon hypersurface. 

The bulk metric is qab- The tangent covector to the horizon is denoted by = (^r, 0^). 
Since is null, it has no canonical normalization. However, we can always match its 
normalization with the choice of coordinates, so that we have ir = 1- Raising the index with 
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the inverse metric, we get the normal vector = g^^is- Being null, it is also tangent to 
the horizon, pointing along its generating lightrays. Therefore, we can restrict the 5d index, 
and write as an intrinsic 4d vector. Once again, this vector does not have a natural 
normalization. It is easy to see that the above conventions fix the following components of 
the inverse metric: 

/- = 0; g^^ = F. (1) 

The bulk metric induces a rich and peculiar geometry on the null horizon. For now, we 
only mention the key concept of surface gravity. Given a choice of scaling for at each 
horizon point x^, it can be shown that the covariant derivative i^Vs^^ is again directed 
along i^. The proportionality coefficient is the surface gravity k: 

i^Vsi^ = /t^^ . (2) 

We stress that k changes under local rescalings of In black hole thermodynamics, 
K, = 2ttT defines the horizon's temperature, and is taken to be non-negative. 



B. Bekenstein-Hawking entropy current 

In this subsection, we briefiy review the concept of a local entropy current in Einstein 
gravity. We start with the Bekenstein-Hawking formula for the entropy of a horizon in 
equihbrium jol, Q : 

S = \a. (3) 

This global area-proportional entropy can always be upgraded into a local current. Indeed, 
any null hypersurface is home to a canonical area current J^rea- This is a vector density 
intrinsic to the horizon and directed along its generating lightrays. As the name suggests, 
its fiux through a spatial slice of the horizon gives the area of the slice. 



•^Area ^e defined as follows 



Start with the spacetime volume density y/—g, and 
dualize it to obtain the volume form ^/^e abode (we consider 5d spacetime as a generic 
example). Now take the puUback of 4 indices into the horizon and contract them with 



g/ivpo-^^l^ where g^^^/"^ is the horizon's Levi-Civita density. The result is a quantity Aa with 
unit 4d density weight, directed along the tangent covector i^- Let us now raise its index 
with the spacetime metric. The resulting quantity A"^ is directed along the horizon's normal. 
For a non-null hypersurface, at this point we could have written A"^ = ^n^^ with the 
unit normal, thus defining the hypersurface volume density ^^7. On a null horizon, however, 
there is no preferred unit normal. On the other hand, the normal A"^ is now tangent to the 
horizon, pointing along its generating lightrays. Thus, we can interpret A"^ directly as an 
intrinsic vector density J^rea horizon. 

By an immediate generalization of ([3]), we can now define a local entropy current as: 

= \JL. ■ (4) 

Does this current satisfy a local increase law d^s'^ > 0? In general, there is no easy answer. 
The quantity 9^ ^Area governed by the twice-projected component Eab^^^^ of the Einstein 
equations. This component is known as the Focusing Equation. When decomposed in terms 
of horizon quantities, it reads: 

ko = Fd^e + + afu) + svrvrr . (5) 

Here, Tab is the bulk energy-momentum tensor; 0"^^^^ is the (non-negative) square of the 
horizon's shear tensor; 9 = d^J'^^.^Jv is the area expansion rate. The "area density" v is the 
coefficient in the collinearity relation J^rea — ^^''j horizon-adapted coordinates, it's given 
hy V = 

We would like the area production rate d^s^ = v6/A to he non- negative. As with the 



global entropy increase AS > 25(|, we see from ([5]) that the answer depends on the null 
energy condition Tab^^^^ > 0. However, this is not enough, since the second-derivative term 
i'^dfj^O can have either sign. This term is related to the fact that an event horizon is defined 
not locally, but by boundary conditions at future infinity which ensure that the horizon 
is eternal. Indeed, these boundary conditions also play a role in the global area-increase 
theorem. 

We conclude that in Einstein gravity, a natural entropy current is always defined, but a 
local entropy increase law d^^s^ > is not always satisfied. In the next subsection, we will 



review an important limit where this law does hold. This (hydrodynamic) limit will occupy 
us for the bulk of the paper. 

C. The hydrodynamic ansatz in gravity 

In non-gravitational systems, an entropy current can be defined in the limit described 
by hydrodynamics. This means that the system is in approximate thermal equilibrium in 
each local region, while the equilibrium parameters slowly vary between one region and the 
next. These local parameters can be encoded by a 4-velocity u^, an entropy density s, 
and optionally also conserved charge densities n"". One can also define the energy density 
e, the pressure p, the temperature T and the chemical potentials The various thermal 
parameters are related by a material-specific equation of state e(s, n"), as well as by universal 
thermodynamic identities. 

A fluid state necessarily requires a separation of scales. On the microscopic side, we 
have the scale determined by the temperature and the underlying dynamics, while on the 
macroscopic side we have the gradients of the equilibrium parameters. In gravity, such a 
separation of scales corresponds to a very large and homogeneous horizon, as compared to 
its surface gravity and the full spacetime curvature. Such horizons are easily constructed as 
planar black holes in AdS. Other examples are also conceivable, e.g. a Rindler horizon in 
Minkowski space. 

To begin with, let us write down a gravitational ansatz for a horizon in global equilibrium, 
satisfying the required separation of scales. We can do so using a constant Lorentzian metric 
h^y and a 4-velocity vector u^, with u^u^ = —1. Here and below, we implicitly raise and 
lower 4d indices with hfj^^, rather than with the gravitational metric g^B- The ansatz for the 
gravitational metric in a neighborhood of the horizon reads: 

9rr = 0; gr^^ = -u^; gf} = /(r)P^^ + k{r)u^Uy , (6) 

where P^^ = h^^ + u^u^ is the projector orthogonal to u^. The coefficient functions /(r) 
and k{r) are constrained by the field equations. Note that with general /(r) and k{r), the 
ansatz holds not just for Einstein gravity, but for any gravitational action which allows 
fluid-like solutions. To describe more general fluids, we may add bulk flelds other than the 
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metric. Indeed, we will consider a bulk Maxwell field in sections llVtlVl The general form 
([6]) of the metric ansatz is again indifferent to these details. The "(0)" superscript on g^y 
stands for "zeroth order" , and anticipates corrections. 

The horizon of the solution ([6]) is located at the value of r where the function k vanishes. 
Without loss of generality, we may take this value to be r = 0. The homogeneity of the 
solution implies that r = is in fact a Killing horizon. Its Killing normal is given by: 

£-^ = (0,w'^); £^ = (1,0^). (7) 

This normalization of is in accordance with the conventions of section III A[ 

The velocity represents the rest frame of the black hole. The nondegenerate metric 
h^y can be interpreted as the metric in which the dual fluid lives. In the AdS/CFT context, 
it corresponds to the metric on the boundary, where the CFT is defined. The horizon value 
/(O) is related to the horizon's area density (recall that A;(0) vanishes): 

Jl.. = V^fiOf/'u^ . (8) 

Thus, the entropy current in Einstein gravity reads: 

= \JU = l^ThfiOf'u'^ ^ V^su^ , (9) 

and we see that /(O) is related to the entropy density s: 

/(O) = (4s)2/3 . (10) 

Similarly, the radial derivative k'{0) is proportional to the horizon's surface gravity, and thus 
to its temperature: 

A;'(0) = -2k = -AttT . (11) 

Unlike the relations (!9|)- (fT0l) between area and entropy, eq. (ITT!) holds also for higher- 
curvature gravitational actions. 

Let us now move on from global equilibrium to one that is local in x'^. We take h^^y, u'^, f 
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and in ([6]) to depend on with very small 9^ derivatives. Note that this is in addition to 
the explicit dependence of / and k on the radial coordinate r, with derivatives which cannot 
be treated as small. Every quantity can now be expanded in powers of the small derivatives 
d^. We will refer to this power as the "order" of the quantity. 

After introducing this dependence on x^, the ansatz (|6]) in general no longer solves the 
field equations. However, it still solves them at leading order. At higher orders, we must 
introduce corrections. By an appropriate choice of coordinates, we can keep the horizon at 
r = 0, and also prevent corrections to Qrr and Qr^- For g^^,, we will consider only first-order 
corrections gj^J. Since we never use the explicit form of gju! in this paper, there is no loss 
of generality involved: g^^J can be understood to contain higher-order corrections as well. 
Finally, we fix the velocity variable u'^ such that ([7]) holds exactly, i.e. always points 
along the horizon's null generators. This constrains the metric correction on the horizon: 

r = : g^^y = . (12) 

In the hydrodynamic limit defined above, the entropy current (jlj) of Einstein gravity does 
have a local non-negative divergence. To see this, let us consider the focusing equation ([5]) 
order by order. Note that the horizon's shear and expansion parameters 6 necessarily 
contain at least one 9^ gradient. Then at first order, eq. ([5]) reads: 

«:(o)^(i) = SvrrilJf^^^ = . (13) 

The last equality follows from our null-energy assumption Tab^^^^ ^ 0, since it's impossible 
to write a non- vanishing positive-definite expression with one 9^ gradient. Thus, at first order 
the entropy production rate d^s^ ~ 6 vanishes. This is always the case with leading-order 
hydrodynamics, known as the "ideal fluid" limit. At second order, the focusing equation 
reads, using 9^^^ = 0: 

^ )2 _^ 87rrig£-^£^ > . (14) 

As promised, at the leading nontrivial order, we see a non-negative local entropy production 
rate d^s'^ ~ 6. The problematic term i^^d^O fell away due to the small-derivative expansion. 
In general, it's not guaranteed that the next-order contributions to (9^s^ will also be non- 
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negative. However, in AdS/CFT, where the metric corrections are controlled by the AdS 
boundary conditions, positivity does hold at the next order; see e.g. |26 |. 



D. Wald entropy 

Wald proposed a definition 12|] for a notion of horizon entropy in general theories of 
gravity. The definition applies to stationary black holes in equilibrium, where the horizon's 
normal (with a suitable choice of normalization) is part of a Killing vector field. The 
entropy is identified with the Noether charge associated with transport along this Killing 
field. The local nature of diffeomorphism invariance then allows us to express this quantity 
as a local integral over a spatial slice of the horizon. 

In the present paper, we consider gravity coupled to gauge fields. We limit ourselves to 
Lagrangians of the form C^qab, Rabcd, ^A, F'ab), with an arbitrary algebraic dependence 
on the Riemann tensor and the field strength. We do not exclude direct dependence on the 
gauge potential Aa, to allow for Chern-Simons terms. In this setup, the rank-2 "Noether 



potential" Q^^ for diffeomorphisms along is given by 



28j: 



= -2£^^^^Vc^D + MdVcC^'''''' , (15) 

where £,^^^^ = dC/ dRAscD is the derivative of the Lagrangian with respect to the Riemann 
tensor. Q^^ is antisymmetric and carries unit density weight. To be precise, in the presence 
of gauge fields, f|T5|) is the Noether potential for "gauge-covariant" diffeomorphisms, where 
together with ^ + ai^ we perform the gauge transformation Aa — )■ Aa — dAidi^As). 
If we considered instead "pure" diffeomorphisms, ( |T5|) would have contained an additional 
gauge-dependent term of the form i'^'Ac ■ dC/dFAB- 

The horizon entropy is now given by an integral of Q^^ over a spatial slice of the horizon: 

S= - I Q^^d^AB = U Q^^'d^AB , (16) 



J T _ 

where k is the surface gravity associated with the horizon normal l^. In Einstein gravity, 
we have C^^'"^ = {y/—g / 327r) {g^'^g^^ — g^^g^'"), and eq. fll6p recovers the Bekenstein- 
Hawking entropy (|3]). We note that when calculating the entropy of stationary horizons, the 
second term in (|T5|) doesn't contribute to the integral [27|, and may be omitted. However, as 



10 



we'll see in section llVt this term is important for the positive divergence of a local entropy 
current in non- stationary solutions. 

III. WALD ENTROPY CURRENT 

A. Definition and possible ambiguities 

We would like to define a local current for the Wald entropy (fT6|) . More precisely, we're 
looking for a vector density on the horizon hypersurface, such that for stationary solutions 
its flux equals the total entropy, while for non-stationary solutions it satisfies d^^s^ > 0. The 
expression f |T6l) for the total entropy suggests an obvious candidate, which is the definition 
we adopt: 

s^' = — g/^'- . (17) 

K 

This formula is written in the horizon-adapted coordinates of section III A[ In coordinate- 
free terms, Q'^'^ is the flux of the 5d antisymmetric density Q^^ through the horizon. The 
result is a vector density on the horizon, with unit 4d density weight. In form language, this 
operation of taking the flux can be described in three steps: 

1. Dualize Q^^ into a 3-form. 

2. Fullback this form into the horizon. 

3. Dualize back within the horizon geometry, resulting in a rank-1 (vector) density. 

Thus, we've defined an intermediate step between the Noether potential Q^^ and the global 
entropy f lT6|) . In f|T6|) . a scalar density was obtained by taking the flux of Q^^ through a 
codimension-2 surface (the spatial slice S). In our definition ( !T7|) . we obtain a vector density 
by taking the flux through the codimension-1 horizon. It's clear that for a stationary solution, 
taking the flux again with respect to the remaining index will bring us back to the scalar 
entropy ( 1T6|1 . 

The definition ( !T7|) can be considered for any gravitational action. Here, we will work out 
its consequences for "algebraic" Lagrangians of the form C{gAB, Rabcd, ^a, F'ab), where 
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qAb jg gjygj^ \yy formula ( !T5|) . Using ( !T5|) and decomposing into radial and tangential 
components, f lT7|) can be written more explicitly as: 

27r 

= _(-2£'^"'^'-(VX- V.4)+44V.£'^""'-) . (18) 

K 

To arrive at this simplified expression, we used the antisymmetry C^^'^^ = _£Abdc ^^^^ 
the relations = and Vj^^i/j = on the horizon. The former relation is due to the horizon- 
adapted coordinates, while the latter follows from Frobenius' theorem £[aVb£c] = 0. We 
left the magnitude of ir in f llSp arbitrary, to exhibit explicitly the dependence on the scaling 
of without entangling it with the choice of coordinates. 

For non-equilibrium solutions, the definition f|T7|) - f|T8|) is ambiguous, and we do not see 
a way to resolve the ambiguity generically. For Einstein gravity however, this can be done, 
as we will describe in section IIII B[ In higher-curvature theories, we must restrict ourselves 
to the hydrodynamic limit. There, at least for actions with an algebraic dependence on the 
Rabcdi the ambiguities can be resolved at the first non-trivial order in small gradients. 

The reason for the ambiguity is that we no longer have a Killing vector field. Therefore, 
we must ask ourselves what exactly is meant by £^ and k in (ITS!) . It is natural to insist 
that remains normal to the horizon, and that k remains the surface gravity associated 
with It remains to determine the local scaling of £^ on the horizon r = 0, as well as its 
behavior at r 7^ 0. The latter is relevant due to the Vr^v contribution to the first term in 
(fTSll . For an unambiguous definition of s^, we must either find a way to fix these freedoms 
or show that they are irrelevant. 

We note that the definition f|T8|) is invariant under global rescalings of i^: the explicit 
factor of is canceled against the implicit factor in k in the denominator. This was of 
course already necessary for an unambiguous definition of the global entropy. Thus, any 
possible ambiguities arise from derivatives of l^. This includes the explicit derivatives in 
the parentheses, as well as the implicit derivative contained in k. 

B. Weak Killing condition 

As a necessary prerequisite, let us now ensure that our definition (fT8|) coincides with the 
Bekenstein current (jl]) for Einstein gravity. In the process, we will fix some of the freedom in 
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the choice of C.^. Having convinced ourselves in the previous subsection that is invariant 
under global rescalings of we use the coordinate freedom to set £r = 1, and look for 
ambiguities only in derivatives of 



For the Einstein-Hilbert Lagrangian, we have: 



32tc 327r 32% 

VX"'"" = , (20) 

= ^v^rr(vx - v,4) . (21) 

Here we encounter a problem. For stationary solutions, the Vjl^ derivative in f l2T]) could 
be related to the Vjy^^ derivative via the Killing condition Va'^b = ~^b^a- To make the 
expression ( 1211) unambiguous, we must find a substitute for this constraint. 



At each point on the horizon, we can fix one degree of freedom for the local scaling of i^, 
plus 5 degrees of freedom for the radial derivative V^^^- Of course, these are not enough to 
force the full Kilhng condition Va^b = —'^b^a- They are enough to force Vr^u = —'^u^r, 
but unfortunately this condition is not covariant: the restriction of a 5d lower index to the 
value r depends on the choice of coordinates outside the horizon. The way out is to force 
the following covariant condition, given by contracting the Killing condition with i^: 

i''{VAiB + VBiA)=0 . (22) 

The /x-component of this condition is in fact an identity: 

^""(V/b + Vb£^) = ^V^(£b£^) + Ki^ = + = 0, (23) 
while the r-component takes the form: 

nVrif. + V/r) = . (24) 

Thus, our condition eats up just one degree of freedom, relating Vr{(^A^^) with the surface 
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gravity. Eq. (12T|) now becomes: 



1 



1 



Ak 



1 

X 



(25) 



In our horizon-adapted coordinates, this coincides with the standard definition 



C. Non-ambiguity at first order 

Returning to higher- curvature theories, let us consider the entropy current (fT8|) in the 
hydrodynamic hmit described in section III CI The current and its ambiguities can now be 
examined order by order in gradients. We can expect better behavior in this hmit than in 
the general case, because is now almost a Killing vector. 

At zeroth order, is a Killing vector, and therefore s'^'-^^ is well-defined. Now, the 
only available 4d vector at zeroth order is u^. Thus, s^'-^^ is always directed along u^, i.e. 
along the horizon's null generators. The only nontrivial content of s^^^^ is therefore its 
magnitude, which is already captured by the scalar area density Q'^^dTiAB- This result is 
in line with the general rule that the leading-order hydrodynamics is completely determined 
by the t/iermodynamics. Another point of view is that at zeroth order, our Wald entropy 
current is simply a rescaled version of the Bekenstein current (jl]), which always points along 
the null generators. 

Thus, the first-order entropy current s^^^^ is where we should start looking for nontrivial 
local information. Here the issue of ambiguities becomes relevant. Once again, we set ir = 1- 
Expanding each piece of (ITSl) in small gradients, the ambiguous contribution at first order 
reads: 

<mbig. = 2vr ^-2£'^-^(o) (^^^'^^'- - ^^^-)) ^ ^ + ^ (^) ^ ^ V.£^'-^'-(°) j . (26) 

It is known [27] that for stationary horizons, the V^jC'^^'^^ term (more generally, the term 
with no factors of V[a^b]) vanishes. The proof is based on the existence of a bifurcation 
surface in the horizon's past Now, the zeroth order in fiuid gradients is locally the 

same as the stationary solution, so for non- stationary solutions we still have = g. 

Therefore, the second term in ( 126|) doesn't contribute. It remains to consider the first term. 
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This requires us to discuss the zeroth-order quantity C^'^'^'^i'^) , The only available rank-2 
tensors at zeroth order are u^u'^ and the transverse projector P^'^ . Thus, C^^'^'^i^) must take 
the form: 

^Mr-r(O) _ ^ (j^p,.. ^ (27) 

The longitudinal coefficient Ci determines the Wald entropy of a stationary solution. The 
coefficient C2 can be shown to vanish, using y^/^Afi^-Co) _ g_ demonstrate this, we 
will need the zeroth-order Christoffel symbols at r = derived from the metric ([6]). The 
non-vanishing Christoffel symbols read: 

KT = -lk'u,, r^ri'^ = ^P!:-, K? = -lu^{f'P., + k'u^u,). (28) 

Note that the Christoffel symbols arise from the radial derivatives of the coefficient functions 
/ and k, which are not small. The vanishing condition for v^£A"~!^^(o) can now be written 
as: 

= VX^"""^^'^ 

_ Y^O) £Arur{0) _|_ pr-{0)£^Ai/r(0) _|_ pi'(0)£/^rAr(0) _|_ pr-{0) £/^r!/A(0) _ pA(0)£^ri/r(0) 

(29) 

pfi(0)^prj/r(0) _j_ pr(0)^/^ri/r(0) _j_ ■pi'(O) ^/^rpr(O) _j_ ■pr(O) ^;iri^r(0) g 

This establishes that C2 = 0. The remaining (ffist) term in fl2B]) can now be written as: 

<mbig. = - 4vrCir (^^r(V.4 - V,4)) ^ ^ . (30) 

Just like in section UlIBl we can now force i^Vriy = —i'^'Vu^r by imposing the weak Killing 
condition (122!) at first order. We are left with: 



(1) 

<mbig. = - 8^C,e^ ( 7^'V.4 j = -SttCiP (^) = . (31) 



We conclude that the Wald entropy current f|T8l) is well-defined at first order. 
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D. Curvature-squared theories 



In the previous subsection, we presented a formal argument for the vanishing of the 
transverse coefficient C2 in fl27|) . This was necessary to demonstrate that our entropy current 
is well-defined at first order. As a consistency check, we will now explicitly show that this is 
the case for curvature-squared Lagrangians. For a given solution, the condition C2 = 0, i.e. 
£ij,rur{o) ^ u^u'^ ^ is linear in the Lagrangian. We can therefore check it for each Lagrangian 
term separately. The Einstein- Hilbert and cosmological terms satisfy the condition, as we've 
seen in section IIII B[ The most general quadratic terms take the form: 

= V—gR^^'^^ Rabcd', ^2 = \/—gR^^RAB', A = y/^R^ . (32) 

Their derivatives with respect to the Riemann read: 

ABCD o / — -jdABCD 



HABCD ^ ]_^^j^ACgBD _ j^AD ^BC _ j^BC gAD ^ rBD ^AC^ (33) 

The relevant components for the question of ambiguities are C^'^^'^; 



-g_ 

£f - = ^gR{g^'g" - g^^ g^ . 



Let us now evaluate these quantities at zeroth order. We will need the volume density, 
the inverse metric and the Riemann tensor derived from the zeroth-order ansatz (161). The 



volume density is ^/--g^ = V—hf^^"^. The inverse metric on the horizon reads: 
while the independent nonvanishing components of R^abcd read: 



(36) 
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Plugging this into we get: 



1 




*-2 



3fk' 



) 



(37) 



2/ 



/^r!^r(0) 
3 



and we find Cf^^'^^'W 



u^u" as expected. 



IV. EINSTEIN-MAXWELL THEORY WITH GAUGE-GRAVITATIONAL CHERN- 
SIMONS TERM 

In section IIII C\ we've seen that our proposed Wald entropy current is well-defined to 
first order in gradients for a class of higher-curvature theories. As the next step, we wish 
to examine the validity of the local Second Law d^^s^ > in a higher- curvature setup. It 
turns out that this task is relatively simple in a 5d theory of gravity coupled to gauge fields, 
with the higher- curvature effect coming from a mixed gauge-gravitational Chern-Simons 
term. The reason for the simplification is that in this setup, the higher-curvature term 
doesn't affect the ideal hydrodynamics, and only shows up in the next (viscous) order in the 
gradient expansion. 

As explained in the Introduction, the hydrodynamics associated with this bulk theory 
is interesting in its own right. The gauge and gauge-gravitational Chern-Simons terms in 
the bulk action correspond to chiral and chiral-gravitational triangle anomalies in the QFT 
underlying the fluid. 

We will flnd that our proposed entropy current for this theory is indeed locally increasing 
at leading order. Furthermore, we will flnd that it coincides with the hydrodynamic entropy 
current in the full fluid/gravity correspondence. 

In this section, we specialize to a fluid living in a flat Minkowski metric, i.e. we set 
h^v = Tj^y. This does not involve a loss of generality. Indeed, it can be shown from index 
symmetries that the axial Chern-Simons terms cannot introduce into our calculations a 
contribution from the curvature tensor TZnuoa associated with the metric huu- That such 
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contributions do not come from other terms in the action was aheady estabhshed in 
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A. Entropy current 

In this subsection, we calculate our proposed entropy current ( IT8l) for the aforementioned 
Einstein-Maxwell theory with Chern-Simons terms, up to first order in gradients. We start 
with the usual Einstein-Maxwell Lagrangian, supplemented with pure gauge and mixed 
gauge-gravitational Chern-Simons terms: 



(38) 



This is the same as eq. (2.6) in p/O], except that A is replaced with —A, and is divided 
by 2 to agree with the normalization conventions of 13l . 

To find the Wald entropy current for this Lagrangian, we write down its derivative with 
respect to the Riemann tensor: 



~ ^e.h. ~^ ^c.s. 

(39) 



r-ABCD _ _ rABCD , rABCD 

On-ABCD 



1 / ^n^[C^D]B 9x Aa (KLMABp CD . KLMCDtj AB\\ 

—99 9 2Aa/i/^ \e Klm + e ^lm ) ) 



167r ' 

The Aa term satisfies the A -H- B and C ^ D antisymmetries and the (AB) o (CD) 
symmetry of the Riemann tensor, but not the Bianchi symmetry. The latter will prove 
immaterial to our definition of the entropy current. We can enforce the Bianchi symmetry 
by redefining C^^'^^ — ). £Abcd _ £[abcd]^ However, we will see in a moment that the extra 
£[abcd] ^gj^jj^ does not contribute to the definition of the entropy current. 

Our entropy current was defined in section IIII Al to be: 

s^^ = — (-2C''"'''iVJr-'^rL)+4:£rV,C''""l • (40) 

Since this definition is linear in C^^'"^, it is enough to show that an antisymmetrized 

term does not, by itself, contribute to the entropy current. This is so, because both L^^'^^ 

and ViyL^'^'^^ vanish when antisymmetrized over the two identical r indices. 
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Due to its linearity in C'^'^'^'^, we can separate the entropy current into two parts: 



c.s. 



(41) 



where ^ is the part of the entropy current coming from the Einstein-Hilbert part of the 
Lagrangian, and s^^ is the part arising from the Chern-Simons terms. The Einstein-Hilbert 
part Sg ^ is given in ([9]). In a flat metric h^u = r]^,^, it simplifies into: 



In general, we define the entropy density s as the norm of the entropy current under the flat 
metric. As we will see, the Chern-Simons contribution s^^ is first-order in gradients and 
transverse to m^. Therefore, it doesn't affect the entropy density at first order. 

We now turn to the higher-curvature contribution s^^ . Using the zeroth-order hydrody- 
namic ansatz (|6]) and the Christoffel symbols found in ( !28|l . one can find by direct calculation 
that R^^}^^ = 0, and therefore CcL'^^^^^ = 0. From the discussion in section IIII C[ it follows 
that there is no zeroth-order contribution Sc.S^ to the entropy current from the Chern-Simons 
term. Our task, therefore, is to calculate the first-order contribution Sc!^}\ 

Before we get started, we must explicitly introduce the gauge potential A'\ into our 
hydrodynamic ansatz. By a choice of gauge, we can force A" = everywhere at all orders. 
With this choice, the zeroth-order gauge potential in the vicinity of the horizon reads: 



The horizon value a''(0,x'^) = —fi°'{x'^) encodes the horizon's chemical potential, while the 
radial derivative 0^(0) = una/ s is proportional to the charge density Ua (cf. [isl). 




Alongside the metric corrections giiu which were discussed in section III (J\ we must also 
consider corrections A^^^-* to the gauge potential fH51) . It turns out, however, that neither 
type of correction affects the entropy current s^^ . To evaluate it, we will need the Riemann 
components: 



(42) 



^a(O) 



a''(r,x^)M^ . 



(43) 




R^iy^ = \u'd,{k'u,) + ^u, 2PP'^d(^u,) + -P'^u'dsf - (/i ^ z/) . (44) 
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From these we calculate: 



r-ixrvrll) _ MQ) f KLM^ir tdW ur , KLMur tdW ^lr\ 

V k' ^ ^ 

Air ZTT 
where the vorticity field is defined as co^ = ^e'^'^P^u^dpUo-- 

We can now calculate the Cc^^^^^\'V u^r — 'V riiyY'^^ piece of (see (HOl) ). At zeroth order, 
the horizon normal £^ = (0,^^) satisfies the full Killing condition Va^b = — Vb^a- Using 
this, we get: 

^M-r-(i)(v^4 - VJuY^^ = 2C>^T^^\V JrT^ = -2£f/'-(i)r::[,°) = -^X^aV , (46) 
where we plugged in the Christoffel symbols from 

We now turn to the second piece of , which is proportional to: 



^ f^-cs. ~^v^c.s. ^ ^ up^c.s. ^ ^ uA^c.s. ^ ^ vp^c.s. ^ ^ uA^c.s. ^ Av^c.b. 
pi nfirvr , -pp, pprvr -pr r<navr , -pr r'prvr , -pr r'prva 

^U^C.S. ~^ Vp C.S. ~^ VCL C.S. Vr C.S. VOL C.S. 



(47) 



As expected, at zeroth order this expression vanishes, since both V"^!^^ and £c.T^^°'' vanish. 
To obtain the first-order expression, we need the quantity: 



£_pupr{0) _ _^j^a{0) ^KLMpv -j^iO) pr _|_ ^KLMprj^{0) 



pu 

LM 



Svr V J 16nf 



(48) 



where we used the Riemann components: 



4^^^ = t^^mK^.' - ^'Pu) - (/i ^ ^) 



(49) 



We note that Cc'^J^'^^^ in (HHjl is antisymmetric in all its Greek indices, so the corresponding 
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terms in (H7|) vanishes. We are left with: 



' ■ ■ ' An Sir 



Plugging f H6|) and fl50|) into (HOjl . we finally get: 

sfP = — T^A^a-^u;^ + — A.aV = -167rTA,/iV , (51) 
K \Ztt Ztc J 

where we used eq. f|TT]) to relate k' to the temperature at zeroth order. Thus, the final result 
for the Wald entropy current reads: 

= s^^ ^ + s^,. = su^" - IQnTXa^i'^uj^ . (52) 



B. The entropy production rate 

As we recall from the beginning of the section, our goal is to demonstrate that the 
proposed Wald entropy current for the theory (!38|) . which we evaluated as (152|) . has a non- 
negative divergence. 

As discussed in section III Bl the divergence of the area current J'Xrea ~ \/~9^^ = Asu^ is 
governed by the focusing equation (the twice-projected component of the Einstein equation 
Eab^^^'^)- In section HTCl we used the detailed form of the focusing equation (|5]) in Einstein 
gravity to prove that the local divergence of the entropy current (jl]) is non-negative at leading 
nontrivial order in the hydrodynamic expansion. 

When we add gauge-gravitational Chern-Simons terms to the Lagrangian, the Einstein 
equations get an extra contribution. This leads to new contributions to the focusing equa- 
tion, which will allow us to prove that the entropy current fl52l) satisfies a local increase 
law. 

We start from the Einstein equation derived by varying fl38|) with respect to the metric: 



^ (Gl + A5f ) = V^g (2F^AcF.^^ - ^F^^f ) 



(53) 



This is the same as eq. (2.16) in [l7|, except that (as mentioned before) F^^ is divided 

21 



by 2, and A is replaced with —A. 

When projected twice along the first row of fl53|) gives us the focusing equation ([5]) 
times a/— with the matter stress tensor Tab = (V^tt) {^F\(jF^^ — (which satisfies 

the null energy condition). 

In appendix [XJ we calculate the contribution of the second line of ( 15^ to the focusing 
equation. This turns out to vanish at zeroth and first order in derivatives. The second-order 
contribution reads: 



new term(2) = 2\jHs {gAcg'''' + S^S^) e^^^^^^V^ {F^klRmnd'') 



(54) 



Thus, at first order the focusing equation ( IT3|) remains unchanged, and we still have 
9^^^ = 0. Using this, the second-order focusing equation reads: 



(55) 



where 9 = d^J'^^^^/ y/—g = 4:d^{su^)/ y/—g. After a short rearrangement, the focusing 
equation becomes: 



> . 



(56) 



We have thus demonstrated that our proposed entropy current (!52|) is indeed locally increas- 



ing. 



V. ANOMALOUS HYDRODYNAMICS FROM HORIZON DYNAMICS 
A. Introduction 

We've seen that a horizon projection of the Einstein equations in a hydrodynamic limit of 
the higher- curvature theory ( l38l) encodes the local increase of the Wald entropy current (l52ll . 
In this section, we will demonstrate that a broader selection of the horizon-projected field 
equations encodes the full equations of charged (anomalous) hydrodynamics. Our proposed 
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Wald entropy current then corresponds precisely to the entropy current as defined in the 
hydrodynamics. 

rn 

The present result extends our previous calculations p^] without the higher-curvature 
Chern-Simons term. As in [l^, the exact correspondence between the horizon field equations 
and hydro equations only holds under certain constraints on the corrections gj^J , A^^^ to the 
fiuid ansatz. In AdS/CFT, these constraints are enforced by the boundary conditions at 
infinity. Apart from assuming the constraints, our horizon-based method does not require 
detailed knowledge of the corrections, because the same hydrodynamics emerges for all valid 
choices of gj},} and A^f!'\ In particular, in an AdS/CFT context, the hydrodynamics which 
we find from the horizon analysis will be the same as the one on the AdS boundary. 

The discussion thus far has been driven by the concept of a Wald entropy current. How- 
ever, the hydrodynamics arising from the theory fl38|) is of interest in a broader context. The 
Xa term in the bulk action generates a peculiar vortical conductivity term ~ T'^u'^, gen- 



erating currents along the vortices in the 4d fiuid. It is conjectured 17|, |l8| that such terms 



are directly related to mixed chiral-gravitational anomalies in the QFT underlying the fiuid. 



In particular, the authors of [iTj] have established this relation in a specific AdS/CFT setup. 
There, the 4d hydrodynamics is embedded in a conformal field theory on the 4d boundary 
of the AdS space. The bulk Chern-Simons coefficients Hatc, Aa correspond respectively to 
chiral and mixed chiral-gravitational anomalies in the CFT. The precise relation is given by 
eqs. (2.9)-(2.14) i n |l7l| (for the conventions in the definition of the anomaly coefficients, see 
eqs. (30)-(32) in [isl): 

^^abc = [tr(TjT,, T,})r - tr(TjT,, TJ)^] (57) 
47r 

A, = --^[tr(T,)^-tr(T,U] , (58) 
where are the generators of the global symmetry group under which the chiral fermions 

n 

in the CFT transform. The authors of [17] then derived the relation between the Xa Chern- 



Simons term and the coefficient of the vortical conductivity term, thus relating the latter 
to the chiral/gravitational anomaly. Our calculation below generalizes the results of jl^ on 
the relation between the bulk Chern-Simons term and the fiuid's vortical conductivity. We 
derive the result for a more general fiuid/gravity ansatz, which allows in particular for an 
arbitrary (not necessarily conformal) equation of state. 
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B. Hydrodynamics from the horizon field equations 



We will now detail the relation between the horizon field equations and hydrodynamics. In 
particular, we wish to recast the field equations into the form of hydrodynamic conservation 
laws. From these one can later extract an entropy current. We begin by writing down the 
Einstein and Maxwell field equations as derived from the Lagrangian ( l38l) : 

V^g [Gl + Ml) = f 2Fl^Ff ^ - ^F^^f ) 

V 2 J (59) 

ds {V—gFa^) = e^^^^^^ ^KabcF^L^MN + -^RrLC^ RmND^^ ■ (60) 



These correspond to eqs. (2.16)-(2.17) in [ij 

The horizon projection of these equations without the \a terms was studied previously in 



13 1 . We briefiy summarize the results. Writing eqs. (I59])-(I60]) symbolically as and M^, 
the relevant components for the hydrodynamics are E^^ and M^, evaluated on the horizon. 
Then, to second order in small gradients, certain linear combinations of these equations 
take the form of hydrodynamic conservation laws for a stress tensor T^^ and a set of charge 
currents J^: 



(61) 



More precisely, the field equations take the form (section V in 13|): 



E^ : -87r(9,T; - J^) = + Ssu^dy/^^ + IQirriau" dy^Af^ (62) 

Ml : AT^d.Jli = . (63) 



The equations depend on the corrections to the fiuid/gravity ansatz through 
the gauge field correction A^*-^^ and the extrinsic curvature quantity c^^"* = —T^fjp'' u^u'' = 
— {1/2) u'^{dyU^ + g'^l^), which appear on the RHS of ( 1^ . The A^j in ( 1^ stands for 
additional terms, which do not fit into the conservation-law template on the LHS, but which 
can be canceled by an appropriate choice of the corrections gjlu^u" and A'^J;^^ (or equivalently, 
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Cu"* and An^^^). With such a choice, we recover the conservation laws in (I6T 



As explained in [13| (see the text surrounding eqs. (V.60)-(V.61)), the correction terms 
cannot themselves produce a contribution to the divergence d^T^. Thus, despite the re- 
maining freedom in choosing the corrections, there is no ambiguity in the stress tensor 
and the charge currents J^, which define the hydrodynamics. 

As calculated in [13] (eqs. (V.21),(V.18),(III.20),(V.63)), the stress tensor and charge 
currents in (l6T]) - (!63|) are given by the constitutive relations: 

s 11^ 16 

T^'^ = euf'u'' + pP"'' n^"" - 2Tfx''aabU^^P''^^dp^ K^b./x^Vw^'^w'') , (64) 

zn T 671 

= nau^ - TaatP^'-'d,^ - -K^kcfi'fi'oo^ , (65) 

In flM|) and fl66l) . the energy density e and the pressure p are derived from the potentials and 
densities T, ^a)S,n°' on the horizon via thermodynamic identities, after assuming an equation 
of state. TT^'" is the shear tensor of the fiuid's velocity, tx^"" = P^'^P'^'d^pU^) - {l/3)P^"'dpUp. 
dab is a conductivity matrix. 

These equations define charged relativistic hydrodynamics, with vanishing bulk viscosity 
and special values for the shear viscosity and for the charge conductivity. The Habc term 
corresponds to chiral anomalies in the underlying QFT. The constitutive relations ( 16^ -( 165|) 
are written in a slightly nonstandard form, since our velocity variable u'^ is neither the energy 
velocity nor the charge velocity of the fiuid. In fact, the horizon normal u'^ corresponds to 
the fiuid's entropy velocity. This will no longer be the case once the higher- curvature 
terms are included. 

Our present goal is to find how the above equivalence to hydrodynamics is modified in the 
presence of the gauge-gravitational Chern-Simons term Aq in the field equations f l59|) -f l60|) . 

The Aa terms on the RHS of the field equations are evaluated in appendices |B] and [Cl 
and read: 

5M: = lQ7x''Ke^^P''dp{Tu,)dp{Tu^) , (68) 
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where the function = p{k'a°' / /)' in the second term in (1671) is unfortunately not express- 
ible in terms of the thermal parameters T, /Xa, s, n". However, we can see from fl62|) that this 
function can be tuned arbitrarily by further adjusting the first-order corrections c}/^ (the 
A^i^^ term can also fulfill this purpose). Thus, by an appropriate choice of corrections, we 
can turn eqs. fl^ - (|^ into: 



(69) 



5M: = lQT,''Ke>'''P^d,{Tu,)dp{Tu„) . (70) 

With this choice of the second term in (169|) . the equivalence ( 16T|) to hydrodynamic conser- 
vation laws is maintained. This is demonstrated in appendix [Dl The constitutive relations 
( |64|) -( 165|) for the stress tensor and charge currents acquire the following extra contributions: 

(71) 

5Jii = -SnX^T^tu". 



For comparison with the literature, it is useful to rewrite the new constitutive relations 
in terms of an energy velocity f ^, instead of the horizon normal u^. The result is: 

Tf"" = (e + p)v''v'' + pri^"" - ^^vr^" , (72) 



TT e + p \6TX 



(73) 



where the (first-order) difference between the velocities and is given by: 

s T e + p \37r / 



These constitutive relations are of the general form given in 2Q| for the hydrodynamics 
of a 4d field theory with chiral anomalies (as derived from entropic constraints). Thus, we 
have demonstrated a full correspondence between our gravitational model and a 4d fluid. 
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The new Ap-proportional terms have the same form as the Pa terms in 20| , which have been 
related in 17|, ll8|] to chiral/gravitational anomahes. 



C. Entropy current from the hydrodynamics 

Given a set of hydrodynamic conservation laws, one should be able to deduce an entropy 
current with positive divergence. For constitutive relations of the form f l72|) - fl73|l . in partic- 



ular with the Xa terms, this has been done in 20|]. In terms of the energy velocity f^, the 
result reads: 

= + ^^a'^'P'^'^d,^ + -^u>' ( A^^^^^^^y _ 27ifian''Xbfi'T) . (76) 
T e + p \3tt / 

Using eq. f l75|) to translate back from the energy velocity to the horizon normal u^^, f l76|) 
simplifies into: 

s'^ = suf" - 167rTAa/i'^w^ . (77) 

This is in precise agreement with eq. ( l52l) . which we derived directly from our definition of 
a generalized Wald entropy current. 

In fact, given the non-negative divergence of the Wald entropy current and the equivalence 
between the gravitational and hydrodynamic equations, it necessarily follows that the two 
entropy currents match. The reason is that, to first order in gradients, a fluid of the type 
discussed in [20] allows just one non- decreasing entropy current to be constructed from the 
hydrodynamic variables. 

VI. SUMMARY AND DISCUSSION 

The local entropy current in Einstein gravity is a vector density directed along the hori- 
zon's generating light rays, whose flux through a spatial slice of the horizon is 1/4 of the area 
of the slice. This entropy current maps via the fluid/gravity correspondence to the hydrody- 
namic entropy current of a field theory. Working in the framework of the fluid / gravity corre- 
spondence, we proposed a definition for a local horizon entropy current in higher-curvature 
gravitational theories. This current is a vector density on the horizon hypersurface, such 
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that for stationary solutions its flux equals the total Wald entropy. We demonstrated that 
our entropy current is well-defined to first order in fluid gradients, for actions with an alge- 
braic dependence on the curvature. It is of interest to carry the analysis further and study 
the properties of the current at higher orders in fluid gradients. 

As an application, we considered the correspondence between 5d Einstein-Maxwell theory 
with gauge and gauge-gravitational Chern-Simons terms and 4d anomalous field theory 
hydrodynamics, which has attracted much interest recently. The chiral and mixed chiral- 
gravitational anomalies in the field theory exhibit new non-dissipative transport coefficients, 
which are responsible for the generation of current along magnetic fields and vortices in the 
fluid. In [20|, the coefficient of the vortical term T'^u'^ in the current came about as a constant 
of integration, and was left arbitrary. It has since been calculated in two regimes: in the 

n ^ 

free-fermions limit [18|, and in the strong-coupling supergravity limit [17]. The vortical 
conductivity was found to be proportional to the chiral/gravitational anomaly, with the 
same proportionality constant in the two regimes. 

In this work, we carried out the analysis of the horizon dynamics along the lines of 



13| . and established fully the fluid/gravity duality between the anomalous quantum field 
theory hydrodynamics and the dynamics of an Einstein-Maxwell black brane with a gauge- 
gravitational Chern-Simons term. In particular, the entropy current as defined in the field 
theory hydrodynamics coincides with our proposed Wald entropy current. The results pro- 
vide further support to the relationship between the vortical conductivity and the mixed 
chiral/gravitational triangle anomaly. In particular, our analysis of the horizon dynamics 
was not restricted to conformal fluids, thus supporting a rather general relation between the 
vortical conductivity and the chiral/gravitational anomaly. Still, it remains to be proven 
that this holds at arbitrary coupling. 
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Appendix A: Chern-Simons contribution to the focusing equation 



In this appendix we calculate the contribution to the focusing equation from the mixed 
gauge-gravitational Chern-Simons term. This contribution comes from the last term on the 
RHS of (|53|1 . twice contracted with the horizon normal: 

(Al) 

Before contracting with u'^, let us write down the full expression for 
{gacSo + Sh9Da) e^^^^^^V E {F^lRmn^^)- It will be needed in appendix O as well. 



(A2) 



— (n A*" -LA''^ \(^CKLMNf) fr^a p DE\.-pC BKLMN j^a p DE\ 

+ e^'^^gsadr {f:,R,^) + e^'^'^gs^d. {F^.R,^^) 

+ r„,,e'^"'^^F;,i?^/^ + 2To.,^ef'''^''F:,R^/' + 2V^,pe^''^^F:^R^r (A3) 

I -pr KXfMU pa p e _|_ rjpr jSKfiu pa p r , rjpr JBkX^i -pa jj r 

_l_ opf" ^/Sk/ii/ pa p e I opr finXfi pa p e 
-t- Zi i^r-'^pi.ya "I" ^J- e/3fc ^ KX-n-^ira ■ 

The necessary Christoffel symbols, Riemann components and field strength components 
are listed in appendix |El 

Since we are working up to second order in derivatives, the first term in (1A3P drops 
immediately, using i^^fj^" = and F^^^^^ = 0. 



When contracted with m", the other terms in the first two lines of (IA3I) also drop, since 
on the horizon we have S'^^m" = 0. 

We are left with the last three lines of ( lA3p . contracted with m". Let's have a look at 
each term separately: 

2u-T%ef^^^^Fj^^^R^^y^ ~ e'^-^'T^u^u- [-u^d,{k'u.) V^u^u^ = 
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a kXiiu pa(l) /pr-(O) p{l) e , pr(l) p(0) e\ _ a kX^u pa(l)pr(0) p(l) e , „,« p(0) e / \ 
z,u, t I -'^Kr ^ ^ rl3 ^ nr I " -^V^^a ~ 

= e^'^^^M, {RiL'u'^Ps) (•••) + r5'^n^(. . . ) = e'^'^'^'^n.n.X. • • ) + = 

= r5^)«^(. . . ) + e^'^^'^r;^ = . 

Summing up, we eventually get: 
2A„m" {g^cSi + 5h9Do) e^^^^-^^Vi, [F^^^Rmn'"') = QAXyTe^'^''d^{fi^u,)d,{Tu,). (A4) 

This is the contribution to the focusing equation from the mixed gauge-gravitational 
Chern-Simons term. 

Appendix B: Chern-Simons contribution to the Gauss law 
In this appendix, we derive the relation (|68|) . i.e: 

e^"'"^i?^.A''i?p.B^ = 327T'e''''P'^d^iTu,)d,iTu^) . (Bl) 

All metric and Riemann components needed for our calculation are listed in appendix |El 
We work up to second order in derivatives. 

We can use the zeroth order Riemann components (eqs. (lElip - (lE13P ) to express 



e^^p" R^yA^ RpaB^ to second order without using the second-order Riemann tensor. We 
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start by writing: 



e i^puA ^pcrB — e'^ ^puAB-tipcr 

= ^'^^"^ {R^LbRT^ + KIabR^J'^ + R'pLbRT^ (B2) 

+ R'^IabR'S'^ + R'^.abR^}'^ + R'XbRT'') + Oid') . 



We see from flET2|) - flET3|) that the contribution R^^^abR^p^^^ vanishes. The R^^^abR?^^^ 



term also vanishes, as follows: 



nupaj^m T3{0)BA _ pvpaj^m d(0) /3« _ /'^' ui^pa p(2) „, 
t ^puAB-'^pa — t J^fiuali-^pa — r ^ ^puajS"''^"- "p 

^ (B3) 



f]^ pupa -d(2) ^ ^13 



In the last equahty, we used the Bianchi identity. For the R^^IarR^p^^^ contribution in (1B2I1 . 
we have: 

piypajDiO) td{2)BA ^ cy avpa td{Q) p{2)ar 
^ ^^puAB^^pcr '^^ ^^pura^^pa 

= 2e--<L + + g'^'-^RT./) (B4) 

The first term vanishes, since R^pvrau'^ = 0. The second term also vanishes, as follows: 

e^^''-<U^"''<i/ = -\fk'e^^'"'uXR%'' = -\rk'e^^P'^u,R^X = , (B5) 

where we again used the Bianchi identity. Putting everything together, we are left with the 
expression: 



^piypcr 



p Bp A _ uupa f r,{0)AB Tj{l) i p(0) p(l)BA 
^puA J^paB — fc I J^ny J^paBA ^pi/AB-^pa 

^ , (B6) 

+ ^R^Xa'^'^^RpI/ + K^BRir^) + oid^) . 

We can now bring the first-order Riemann tensor to the standard index placement Rf^uA^ ■ 
The indices of R^^Ia^ can be raised and lowered with the zeroth-order metric. When raising 
or lowering R^^Ia^ i we get expressions of the form g^^^R'^^'' + g'^^^'R'^^'' (with indices sup- 



31 



pressed). Whenever an R^^^ factor is multiplied by another it can be dropped, since 

?(0) Bp(0) i 



g^i/p(T^(o)^B^(o) D _ g antisymmetric product of two Uu's. In this way, the first 



two terms in ( 1B6P can be written as 26^"'"^ R^^^y^^ R^pJ^"^ . For the last term, we have: 
Using gfW^^^^ = and QaW^^^^ + 9aW^^^^ = 0, this can be written as: 

uup(7tdW o{1)BA _ nupa f B d(1) A p(0) ^(1)SC d(1) A d(0)BA^(1) p(l) c\ mo\ 

<^ ^puAB^pa — ^ y^puA ^paB ^puAB iJ ^pcrC ^pv SBC^paA J ■ ) 

Putting things together again, we end up with: 

_ ^(I)^C'r(I) ^_ P(0)B^^(1) d(1) C I p(l) Bp(l) A 

^pvAB y ^pcrC ^pu yBC^paA ^puA ^po 



"-puAB y ^"-paC ^"-pu yBC^'-paA ^ puA paB 

pupa (nviO) 13 oil) A, p(0) Jl)aA r 
^ y'^^pfA ^paP + ^puraS ^paA 

_ f?iO) °_ /?(0)/3a„(l) d(1) C , p(l) B p(l) A 

_ p(0)/3a (l)p(l) A , p(l) Bp(l) A 



(B9) 



The Riemann components R^^Ia^ ^^n be found in eqs. (1E21I) - (1E27P of appendix [El 

We are now ready to evaluate the expression (]B9p . For the R^^Ia^R^p1i3^ term, the relevant 
components of R^p}p^ are R^p^i^^ and R^^^i^^u^ . Their product with ^'^"^'^ R^^^a^^ vanishes, 
due to the antisymmetrized product of two m^'s. For the second term in (IB9p . we have: 

e''''''R^X{9^'^''^K'^A' + y^'^'^R^^^^^ , (BIO) 

where h^i, is the fluid (flat) metric. This term again vanishes due to an antisymmetrized 
product of two "U^'s. For the third term in (\B9\i . we have: 

^..p.^(o)/3a^(i)^a)^A ^ e-^p'^^oj/.^^w^a)^. ^ ^e'^-^'Val'^lR^l' ■ (bh) 

Again the relevant components are R^pJa'^u'^P^, whose product with e^^'^^u^, vanishes. We 
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are left with the product ^^^'^'^ R^^Ia^ ^^^Ib^ ■ This decomposes as: 

nupa -Bp(l) A _ uupa ( p(l) r p{l) r . n Ty{l) a r . p(l) Z?,,",, „ p(l) ^.,1.,^ 

(B12) 

The e^'^'P'' R^p}/ B^-p}/ and e'^'^'^^i^il^a^M^M/ji^ii-VM^w^ terms contribute e^"'P''df,{k'u^,)dp{k'u„) 
each. The other terms don't contribute, as we now show: 



We finally get: 

^^'^'"^ R^pIa^ R^puB^ ~ e^'^'"^ {R^pu/R^pu/ + R^pua^'^°''^l3R^pJ-y^u"'Us) 

e^^P'-R.^A'^RpuB^ = e^-'^'-R^^A^'K^^ + 0{d') 
Eq. (IBip is therefore demonstrated. 



(B13) 



(B14) 



= 6'^'^^'^ + =0 + = 

e'^'^''-i?«/t."Pj/?« = e'^^^'^i^W/^-p; (t^p (...) + ^P7 (•••)) 

= eP-'P'^ {u,u, (...) + P«/«"P,, (...)) 

= + k'e^-'P'-u, QjPupU^dxf + P,^5(An,)^ (...) 
= + k'e'^'P'^u^ (^d^.up) + iw^M^^AM.^ (...) = + = 

^,.p.^(i)^^pap7^aM = e'^^P^ (p^^ (...) + p^^ (... )) (p,^ (...) + p; (... )) 

= 6'^'^^'^ (p,,(...) + p,„p,,(...) + p^^p; (...)) 

= + e^-''- (/'2p;P^^P;P^^9(An^)9(,M5) (B15) 

+ p^^p; (p..w^9a/ + /p;pX^a)) (...)) 

= + e^^^-^P; {P^yOJ + fP^P^d^^u^)) (...) = . 



(B16) 



(B17) 
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Appendix C: Chern-Simons contribution to the Gauss-Codazzi equation 



In this appendix, we derive eq. (l67|l . i.e.: 



\9t^C9 +OcO^)e y E {i" KL^MND ) 



(CI) 



where Q'^ = (a"(/F - f'k') + a'^fk') = fik'a^/f)'. 

In appendix El eq. (IA3p . we have already written the full expression for the LHS of flUip . 

Let us now use the Christoffel symbols, Riemann components and field strength compo- 
nents from appendix[E]to evaluate each term in ( lASp separately (this time without projecting 
their vector index). 

Up to second order in derivatives, we have: 



Uad^iu^a'^-k') + a'°-k' yu^d^uJa - -Uad^u' 

a fk" f'k' 



2 4/ 



-AfP^X^^^d^{a'^ux)d^(^-^ 
+ {2k" f - f'k') a'* [2uJa,d,u' + 2u'd,Uo, - 



f'k' 

ia'^k')' - a^'^-j- 



+ 2fcu. 



'^"^ oo i-A a"'k'f' a-^ r 
U Oaf ] U dc^f 



f 



P 



2k'f' 



+ e^"^'^a.(a'^MA)P^« ( -^dj + f'd.k' + 2k'd,r 



e^^'-'gasd, {f:^R,J') = -2fuad, (""'^T^) + 2aV'fc'w^5e«a - a'' f'k'{2u'd,Ua - u^d, 
Ta^re^'^'F^^R,/' = k'uae^'^''d^ia''ux)d,ik'u,) - 2a'' f'k' (^u'd^Ua + ^UaU^f^ 
2Tae^e^''^''F:,R,/' = 2a"'k'uaU,dJ 
2V„,Be^^^^F:^R,r = 2o^ ( k" - ^] u^u^dj 



V 
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a(0,l) 



— P 



2j.,^^^,^,^,^^^^^, = 



where we have used the following definitions: 



«/3 + afsU^ - UpU^Qsl^u" 



acceleration vector : 



vorticity vector 



(C2) 
(C3) 



as well as some geometrical identities (explained in section V of [13]): 



(C4) 



Summing all the terms, we get: 



2k'uae''^^''d,{a''ux)d^{k'u,) 



(C5) 
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where we defined Q°'{x^) = {a"-{fk" — f'k') + a'^-fk') = p{k'a"- / f)' , and used the definitions 
flE2p of the thermal parameters. We also used the the first-order focusing equation d^{su^) = 
(entropy conservation) for the last equality. Eq. (ICip is therefore demonstrated. 



Appendix D: Chern-Simons corrections to the hydrodynamic constitutive relations 



In this appendix, we derive the contributions f[7T|) to the charge currents and the stress 
tensor from the Xa terms in the Einstein-Maxwell equations (!70|) - (!69|) . 

Using the equivalence fl6T]) to relate the new contributions to conservation laws, one gets: 



-i^5E; + i^/iX5M: a.(5T;)=0 ^^^^ 



An 



6m: d,i6Ji:) = o. 



Then using the explicit form of the new contributions ( !70l) -(!69l) (remember those were on 
the RHS of the field equations), one gets: 

-8TrXaS^2Tu^d,{fi''Tuj^) - Asu^di^ (^'^'^l) + f^'u^T^^l^ ^ -^.T; = ^^^^ 



where we have used 



e'''""'d,iTu,)d,iTu^) = 2d,{T^uj^) (D3) 

e^''P''d^{li''u,)d,{Tu,) = 2d,{^''Tun . (D4) 

For the stress-energy tensor, some more work is required. Let us develop it a bit: 

- 8nXaS^2TuMl^"Tu'') + fi'^uMT^^") ' '^■^u''d[^ (^^^lu,^ | = (D5) 

- SnXalsd^ifi^T^u") - 2d,{Tu^)fi''Tu;'' - d,{fi^u^)T^u'' (D6) 



- Au'di^ifi^'T^u,^) + 2suVT2u;^9, | . (D7) 

Using the first-order focusing equation d^{su^) (see text above eq. (|T^ ) and remembering 
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we are only interested in second order in derivatives, this becomes: 
Using the hydrodynamic identity: 

1 1 

(^"diuUf,] = -Uf^oj^ai, = -u^duUj" , (D8) 
where = u^dj^u^ is the acceleration vector, only the first term survives. We finally get: 

ending up with: 

ST^"^ = -327rA„^"T\('^cj'^) (D9) 
5J^ = -MXaT'^uj" (DIO) 

We have thus demonstrated eqs. ( ITTi) . 



Appendix E: List of ChristofFel symbols, Riemann components and field strength 
components 

We now summarize the Christoffel symbols, Riemann components and field strength 
components which we need for the calculations in the previous appendices. 

As mentioned in sections III CI and IIV Al we use the following hydrodynamic ansatz for 
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the bulk fields (metric and gauge potential) in the neighborhood of the horizon: 



(0) 

(El) 



grr = 0; gr^, = -u^] gf} = f{r, x^)P^y + fc(r, x^)u^u 



with first order corrections to the metric gju! satisfying g'^^Ju" = on the horizon, and first 
order corrections A^*-^-* to the gauge potential. 

On the horizon g^^,^ becomes degenerate, i.e. k{0) = 0. The other functions and their 
radial derivatives are related (at zeroth order) to thermal parameters as follows (cf. Ill C\ 

/(O) = (4s)'/3; A;'(0) = -AnT; a"(0) = -/i"; a"^(0) = nnys , (E2) 

where primes stand for derivatives with respect to the r coordinate. 

The inverse metric in the neighborhood of the horizon at zeroth order immediately follows: 

The corrections to the inverse metric are found as: 

gWAB ^ ^giO)ACgWgmDB ^ ^^^^ 



This gives: 



where "r = 0" indicates the horizon value, and h^^ is the 4d fluid (flat) metric. 
The conventions for the Christoffel symbols and the Riemann are as follows: 

^ABC = d{Bgc)A - ^dAgBc (E6) 
38 



^BC = g^^^DBC 

Rabc^ = dsT^c + ^BE^AC - {A'H- B) . 



The non- vanishing lower Christoffel symbols at zeroth order read: 



The non- vanishing upper Christoffel symbols at zeroth order are: 
1 /' 1 



The Riemann tensor components of type R^va^ read: 



The nonzero components of R^^Iab are: 



R^^l,^-\fk'l^P.,-{lJi^v) 



The Riemann tensor components of type R^^"^^ are: 



Some more Riemann components at zeroth order which we will need: 



d{0) r ^ „(0)r/3 (0) d(0) P _ i^p 



ap ■ 
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The first-order lower Christoffel symbols are constructed from the derivatives d^g^^^^ and 
^rQAB- The nonzero components are: 

1 1 

^% = - i^P^P^nf + Pn{vdp)f - -u^Upd^k + u^U(„dp)k (E15) 
+ if + k) {-dp,U(^yUp) + dpU^i^u,,) + d^U(nUp)) . 

The first-order upper Christoffel symbols are found as: 

r^? = 5^^^^^rgic + ^^°^^^rgic- (ei6) 

The nonzero components read: 

= + - (jU^g'^JP^ ^l{a, + u^g'S)) (E17) 

rM^) = i + \p^'g'^^ - (^P^'glSPu) (E18) 

^ ^ ^ ^ (E19) 

- U(,j.dv)k - fPPP^d^pU^) - ka^f,u^) --P^^u^dpf - fP^P^d^pU^) 

n^p^ = J (-lP^^pP''%f + P(M - lu.UpP^'^d^k^ 

+ + {-P^'^daU^^Up^ + U^^dp^U^) - (d^^Up) + i^^)^ 

+ ^^'alSP'' ifP^P + k'u.up) ^ i (^Pi;dp)f - ^P^pP'-'-dj'^ 

+ U^^dp^ul' - dl'u^^Up) - U'^ {p^PXaUx) + \g'^p^ , 

where 

a" = u^d^u" 

is the acceleration vector. 

First-order Riemann components whose r dependence is needed: 

4^' = Im^'uu) + Yf {kd[uUtA - \u,d,k^ -ifx^u)^ Id^ik'u,) - {f^ ^ i^) (E21) 



(E20) 
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+ 



2 

kf_ 

V 



kf f k\ kk! 



r=0 



1,, 



> - -k'u^ [ -Pr.pU'^dJ + fP,^P;d^^U„) 



^'J'" = lu^u^d^k' + (k' - \^u^d,u, + i (^1 - I ) u^d.u" 



f 



k 



- ^d^kiuf'u. + pp) + Aa^/' pp 



irPH ( k' + ^]dj 



f 



r=0 1 k' 

^u^dp,{k'ui,) + —Ui, 



2Pf"'d^,u^) + jP^u'dsf 



(E22) 



(E23) 



For the other needed Riemann components (whose r dependence is not needed) , we have: 



- k' (^2P'''d[,u,] + Po^g'^l^ - ^P^^^W^ } 



-^i^y/fP^pU^d,^ + UpdM + k' {d^^Up^ - ^ I + fPid{xU.) 
+ ^ii'^ f - fP^pP^dJ^ 



f 



+ Ifp^ (JjP^p^'d,f + p;p,^a(.«,)) - (/^ ^ 



lira 



+ -[fP^a + k'u^Ua]u'g'^l\'' 



(E24) 



(E25) 



- ^Pa(.9^^^.) - \ VfPa^u'd, (^) + J {u.u'gS + P^O (E26) 
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Sometimes it's useful to decompose the last two indices of R^M^p"^ along and orthogonally to 



R'^pl'^PpU. ^ \up (^^/fP.pu'd,^ + fP^d^xUu) + k' {P^dy^u,^ - \gf^Pi 
R'pL^PpP^ ^ If [Pp (^jP.pu'dxf + P^P^d^.u,^^ - P^.P^P'^^d^^u,-, 
As for the field strength tensor, all we will need is: 

^a(O) ^ ^a(l) ^ Q. ^a{0) ^ ^a^^^ ^p-^^^ 



(E27) 



(E28) 
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